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We  present  an  extension  of  the  rigorous  coupled-wave  analysis  technique  to  analyze  second-harmonic  genera¬ 
tion  (SHG)  in  periodic  optical  nanostructures  in  the  undepleted-pump  approximation.  We  apply  this  method 
to  analyze  SHG  in  two  example  nanostructures  for  which  we  predict  enhanced  nonlinearity  due  to  transverse 
near-field  localization  of  the  fundamental  optical  field  in  the  nonlinear  material.  First,  we  examine  a  periodic 
nanostructure  that  yields  up  to  twice  the  transmitted  SHG  intensity  output  compared  with  the  bulk  nonlinear 
material  but  only  for  small  nanostructure  depths  because  of  mismatch  of  the  fundamental  and  second- 
harmonic  mode  phase  velocities.  Second,  we  develop  and  analyze  a  modified  nanostructure  and  find  that  this 
nanostructure  concurrently  achieves  transverse  localization  and  phase  matching  for  SHG.  In  principle,  this 
permits  an  arbitrary  coherent  interaction  length,  and  for  several  specific  nanostructure  depths  we  predict  a 
transmitted  SHG  intensity  output  more  than  two  orders  of  magnitude  greater  than  that  of  the  bulk  material. 
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1.  INTRODUCTION 

Recent  developments  in  microfabrication  techniques  and 
nanoscale  photonic  technologies  suggest  that  integrated 
optical  or  photonic  systems  based  on  these  technologies 
may  be  realized  in  the  near  future.  In  particular,  any 
high-level  integrated  photonic  system  will  almost  cer¬ 
tainly  require  the  inclusion  of  nonlinear  optical  processes 
(e.g.,  to  achieve  wavelength  conversion,  switching,  signal 
amplification,  or  system  reconfigurability).  However, 
most  standard  nonlinear  optical  materials  are  not  readily 
compatible  with  standard  microfabrication  materials  and 
techniques.  Fortunately,  nanostructured  material  and 
device  technologies,  including  artificial  dielectric 
nanostructures,1,2  form-birefringent  nanostructures,3,4 
and  photonic  crystals, 5-7  may  provide  a  means  to  over¬ 
come  these  incompatibilities.  Since  optical  nanostruc¬ 
tures  permit  the  engineering  of  their  effective  optical 
properties,  they  facilitate  realization  of  the  desired  linear 
and  nonlinear  optical  functionalities  with  use  of  compat¬ 
ible  materials  and  processes. 

Photonic  nanostructures  not  only  facilitate  the  integra¬ 
tion  of  various  optical  devices  but  may  also  permit  im¬ 
provement  or  enhancement  of  their  optical  functionality. 
We  have  previously  investigated  transverse  localization  of 
the  optical  field  inside  a  periodic  nanostructure  and  its 
potential  application  to  the  enhancement  of  nonlinear  op¬ 
tical  phenomena  such  as  second-harmonic  generation 
(SHG).8  We  predicted  that  the  electromagnetic  fields  in¬ 
side  the  nanostructure  would  experience  strong  trans¬ 
verse  localization,  significantly  elevating  the  field  ampli- 
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tude  in  the  nonlinear  material  and  thus  enhancing  the 
SHG  output.  Similar  enhancement  of  nonlinear  optical 
effects  though  near-field  localization  in  a  variety  of  sub¬ 
wavelength  optical  structures  has  been  predicted  and  ob¬ 
served  by  a  number  of  investigators.9'13  Further  nonlin¬ 
earity  enhancement  may  be  realized  by  achieving  phase 
matching  between  the  fundamental  and  second-harmonic 
(SH)  fields  in  a  nanostructure.14'19  To  achieve  the  larg¬ 
est  SHG  output,  we  investigate  a  nanostructure  that  em¬ 
ploys  both  enhancement  mechanisms:  near-field  local¬ 
ization  and  improved  phase  matching. 

In  order  to  accurately  design  and  analyze  nonlinear  op¬ 
tical  processes  in  subwavelength  nanostructures,  it  is 
necessary  to  develop  appropriate  rigorous  modeling  tools. 
The  rigorous  coupled-wave  analysis20  (RCWA)  technique 
is  a  well-established  tool  for  analyzing  linear  electromag¬ 
netic  wave  interactions  with  subwavelength  periodic 
structures  in  both  the  near-  and  the  far-field  regimes.  In 
this  paper  we  extend  the  RCWA  method  by  using  a  per- 
turbational  expansion  of  the  nonlinear  polarization  term 
of  Maxwell’s  equations  to  analyze  the  nonlinear  SHG  pro¬ 
cess  in  the  undepleted-pump  approximation.  The  result¬ 
ing  electromagnetic  analysis  tool  permits  an  accurate  in¬ 
vestigation  of  the  SHG  enhancement  in  periodic  optical 
nanostructures.  In  Section  2  we  describe  the  extension 
of  RCWA  to  model  SHG  in  the  undepleted-pump  approxi¬ 
mation.  In  Section  3  we  apply  this  tool  to  study  the  SHG 
enhancement  resulting  from  two  nanostructures:  the 
original  design  of  Ref.  8  and  a  modified  structure  opti¬ 
mized  for  phase  matching.  We  present  a  summary  and 
conclusions  in  Section  4. 
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2.  ELECTROMAGNETIC  ANALYSIS 
METHOD 

A.  Frequency-Separable  Formulation  of  Maxwell’s 
Equations 

To  analyze  SHG  in  the  undepleted-pump  approximation 
by  using  RCWA,  we  assume  that  electromagnetic  fields 
exist  only  at  the  fundamental  and  second-harmonic  fre¬ 
quencies,  yielding  total  fields 

E(r,  t)  =  EJf)exp(jat)  +  E2ol(r)exp(j2wt),  (1) 

H(r,  t)  =  H^expijwt)  +  H2lo{r)exp(j2ut),  (2) 

where  E w{r)  and  E2co(r)  are  time-harmonic  representa¬ 
tions  of  the  electric  fields  at  the  fundamental  and  second- 
harmonic  frequencies,  respectively,  and  H (J  f  )  and  H2lo(r) 
are  time-harmonic  representations  of  the  magnetic  fields 
at  the  fundamental  and  second-harmonic  frequencies,  re¬ 
spectively.  We  assume  p,  =  p0  for  simplicity.  Since  we 
wish  to  investigate  nonlinear  phenomena  in  nanostruc¬ 
tures,  we  must  allow  the  electric  displacement  to  include 
both  linear  and  nonlinear  effects: 

D(r,  t )  =  eE(r,  t )  +  PNL(r,  O-  (3) 

However,  for  simplicity  we  assume  that  the  nonlinear  po¬ 
larization  term  contains  only  the  sum-frequency  genera¬ 
tion  of  the  second  harmonic  from  the  fundamental  fre¬ 
quency 

x(2)  - 

TNL(fi,  t)  =  —  Elll(r)exp(ju)t)Elll(f)exp(j^t)  (4) 

and  ignore  all  other  y':2:  phenomena  such  as  difference- 
frequency  generation,  optical  rectification,  or  cascaded 
nonlinear  effects,  as  is  consistent  with  the  undepleted- 
pump  approximation.  Substituting  into  Maxwell’s  equa¬ 
tions,  we  find  that 

V  X  [EJ^exp^jcot)  +  E2lo(r)exp(j2tot)] 

d  „ 

=  -fi  —  [HJf)exp(jcot)  +  H2ol(r)exp(j2cot)]  (5) 

ot 

V  x  [HM(r)exp(j(ot)  +  H2ul(f)exp{j2u>t )] 

d 

=  ~  [ewEco{f)exp(j(ot)  +  £2(I>E2ol(f)exp(j2iot) 

ot 

+  TNL(r,  t )],  (6) 

where  em  and  e2w  are  the  dielectric  constants  of  the  mate¬ 
rial  at  the  fundamental  and  second-harmonic  frequencies, 
respectively,  which  in  general  will  be  different  owing  to 
material  dispersion.  Equations  (5)  and  (6)  can  be  decou¬ 
pled  into  two  sets  of  coupled  differential  equations  at  the 
fundamental  and  second-harmonic  frequencies.  At  the 
fundamental  frequency  we  obtain  the  standard  homoge¬ 
neous  form  of  Maxwell’s  equations: 

V  X  Ew{r)  =  -jwpHm(r)  (7) 

v  X  HJr)  =  ;»ew£u(f).  (8) 

The  propagation  of  the  fundamental-frequency  field 
through  the  structure  is  determined  by  using  the  stan¬ 
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dard  RCWA  approach  and  is  summarized  in  Appendix  A  to 
introduce  the  notation  used  in  the  remainder  of  this  sec¬ 
tion  as  well  as  to  provide  a  self-contained  description  of 
the  SHG  RCWA  method.  At  the  second-harmonic  fre¬ 
quency,  however,  we  obtain  an  inhomogeneous  system  of 
equations  that  are  due  to  the  nonlinear  polarization  term: 

V  X  E2co(r)  =  -j2wpH2ol{r),  (9) 

V  X  H2m{r)  =  j2ioe2loE2ltl(r) 

+  jo>X(2)EJf)EJr).  (10) 

This  inhomogeneous  system  of  equations  can  be  treated 
with  standard  techniques  and  is  described  in  detail  in  the 
following  sections. 

It  is  important  to  note  that  useful  information  about 
the  eigenmodes  of  the  periodic  structure  can  be  readily 
obtained  from  this  formulation  of  the  RCWA  method. 
The  electromagnetic  fields  inside  the  structure  are  repre¬ 
sented  as  a  superposition  of  the  eigenmodes,  as  shown  in 
Eq.  (A8).  Each  column  of  the  eigenvector  matrix  Wm  de¬ 
scribes  one  of  the  eigenmodes  of  the  grating,  with  the  re¬ 
spective  eigenvalue  corresponding  to  the  propagation  con¬ 
stant  of  that  mode  and  the  coefficient  corresponding  to 
the  complex  amplitude  of  that  mode.  This  information 
can  be  used  to  identify  the  energy-carrying  modes  of  the 
grating  structure  and  to  compare  the  propagation  con¬ 
stant  and  the  transverse  profile  of  each  of  those  modes. 

B.  Nonlinear  Polarization 

In  the  standard  RCWA  formulation,  owing  to  the  periodic 
nature  of  the  grating  a  space-harmonic  expansion  of  the 
dielectric  constant  function  of  the  grating  is  performed. 
In  a  similar  fashion,  we  can  also  expand  the  nonlinear  co¬ 
efficient  x‘l  )  using  a  space  harmonic  expansion: 

X{2\x)  =  2  Xk  exp(jkKgx).  (11) 

k 

Using  the  expansions  shown  in  Eqs.  (A2)  and  (11),  we  can 
write  the  undepleted-pump  nonlinear  polarization  term 
as 

P NL*r)  “  2  S  2  XkS<o,iSv,h 

i  h  k 

x  exp{-j[(dWii  +  dW|/l)  •  f  -  kKgx]}. 

(12) 

Expanding  the  argument  of  the  exponential  function,  we 
find 

(<? u,,i  +  •  r  -  kKgx 

=  [S&inc,^*  ~  (i  +  h  +  k)Kg]x  +  2kiwmyy.  (13) 

Since  the  fundamental  field  and  the  nonlinear  coefficient 
are  all  represented  by  harmonic  expansions  based  on  the 
period  of  the  grating,  it  is  possible  to  write  the  nonlinear 
polarization  term  as  a  single  space-harmonic  expansion, 

■Pnl,u(u  t)  =  2  P,,v(z)exp{-j\(2kinCtx  ~  iKg)x 

i 

+  2  kincyy]},  (14) 

where 
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Pi,u(z)  \-Xvxx,i—h—k^*<o,h,x^*a>,k,x 

h  k 

Xvxy,i~h—k^io,h,x^(o,k,y 
Xvxz,i—h—k^io,h,x^(o,k,z 
Xvyx,i—h—k^a),h,y^(o,k,x 
Xvyy,i—h—k^(o,h,y^(o,k,y 
Xvyz,i~ h~ io,h,y^ (o,k,z 
Xvzx,i—  h—  io,h,z^ oj,k,x 
Xvzy,i—h—k^io,k,z^io,k,y 
Xvzz,i—h,—k^(o,k,z^io,k,z\  (15) 

and  u  g  {x,  y ,  z}. 

In  analyzing  the  fields  at  the  second-harmonic  fre¬ 
quency,  we  perforin  a  space-harmonic  expansion  of  the 
fields  similar  to  that  of  Eqs.  (A2)  and  (A3): 

E2 Sr)  =  2  {S2u,,i,x(z)x  +  S2oi,i,y(z)y 

i 

+  •S2w,i>3(2)2}exp(-J<?2,U;;  •  r)  (16) 

H2ol(r)  =  A/ — 2  {U2a,i,x(z)x  +  u2«,,i,y(z)y 
v  Vo  i 

+  U2a,'i,z(z)z}exp(-ja2l0!i  ■  f)  (17) 

To  find  the  transverse  wave  vector,  we  consider  that  the 
second-harmonic  field  is  generated  by  the  fundamental 
field,  and  thus  the  transverse  components  of  the  second 
harmonic  wave  vector  must  be  directly  related  to  those  at 
the  fundamental  frequency: 

&2w,i  =  (2&in<vc  -  iKg)x  +  2 kincjy.  (18) 

To  match  the  range  of  the  space-harmonic  expansion,  the 
second-harmonic  field  expansion  must  have  2p  —  1 
=  4a  +  1  components  indexed  by  i  e  {-2a,...,  -1,0, 
l,...,2a}. 

C.  Solution  of  the  Second-Harmonic-Frequency  Fields 

Using  the  inhomogeneous  equations  for  the  fields  at  the 
second-harmonic  frequency,  and  inserting  the  space- 
harmonic  expansions  for  the  electric  and  magnetic  fields 
shown  in  Eqs.  (16)  and  (17),  respectively,  as  well  as  the 
expansion  for  the  nonlinear  polarization  shown  in  Eq. 
(15),  we  obtain  a  coupled  system  of  equations  in  matrix 
form: 


j  S2m,x(z)\ 

I  S2oljX(z)\ 

d 

S2a,y(z) 

—  A-2  (j) 

S2o,,y(z) 

dz 

U2(0fx{z ) 

^2  w,x(z  ) 

\  ^2^)1 

\U2w,y(z)l 

™2  io,x  e2a),zzJr  z* 

^2a),ye2oj,zzPz(z) 

Py(z')  ^2io,yz^2at,zz±  z[ 

\  Px(z)  T  ^2ai,xz^2(o,zzP  ziz) 


zPz(Z) 

Z) 

zPz(z) 


where  the  matrix  A2w  represents  the  coupling  between 
modes  that  are  due  to  the  properties  of  the  grating  and 
Maxwell’s  equations  and  is  essentially  the  same  as  in  the 
homogeneous  fundamental-frequency  solution.  This  in¬ 
homogeneous  system  of  equations  can  then  be  solved  by 
using  standard  techniques.  Defining 


Sjm,x(z)\ 

^  2(ii, y(z^ 

(o,y(z)  I 


=  W  2uQ(z), 


(20) 


where  W2(u  is  the  eigenvector  matrix  obtained  from  ma¬ 
trix  A2u)  [equivalent  to  matrix  W,„  shown  in  Eq.  (A8)  but 
for  the  second-harmonic  frequency),  we  can  write  the  sys¬ 
tem  of  partial  differential  equations  in  Eq.  (19)  as 


-Q(z)  =  W2w1A2wW2wQ(z)  +  W2({Y(z), 
dz 


(21) 


where 


Y(z)  = 


Z2(i),x^2o),zzJr  z 


k 


2u),yt2(o,zzjrz 


Pz(z) 

Pz(z) 


(22) 


Py(z)  ^2(o,yz^2ai,zzP  z^z) 

Px(z)  ^2oj,xz^2co,zzP Z^z)  I 

The  partial  differential  equation  shown  in  Eq.  (21)  con¬ 
sists  of  a  homogeneous  and  an  inhomogeneous  part,  so  we 
explicitly  separate  the  homogeneous  and  the  inhomoge¬ 
neous  parts  of  the  solution: 

Q(z)  =  Qih\z)  +  Qm(z).  (23) 

The  homogeneous  part  of  the  solution  is  the  same  as  in 
the  fundamental-frequency  case: 

Qlh)(z)  =  c2(Ji,  exp[\2(J;!(z  “  2<m)]»  (24) 

where  c2tu ,  is  the  amplitude  coefficient  of  the  ith  eigen¬ 
mode  to  the  determined  by  the  boundary  conditions,  \2w  ; 
is  the  eigenvalue  corresponding  to  the  ith  eigenmode,  and 
2o,i  is  defined  by 

\d,  0  =£  i  <  2p 

Zo,i  =  ( 0,  2p^i<4p  (25) 

The  renormalization  shown  in  Eq.  (25)  is  performed  to  en¬ 
hance  numerical  stability  (see  Ref.  21  for  a  more  detailed 
explanation).  The  inhomogeneous  part  of  the  solution 
produces  the  inhomogeneous  term  in  the  partial  differen¬ 
tial  equation  of  Eq.  (21)  and  is  given  by 


exp(\; 


„,iz)  P 

Jd 


exp(-\2  ;z') 


Q!:ih)(z)  = 


(19) 


x  S  W2ii0-y/z')dz',  0  «  i  <  2p, 

j 

exp(X2w,;z)  exp(-\2wii2') 

J  0 

X  2  W2i^(z')dz',  2p  «  i  <  4p, 

\  J 

(26) 

with  use  of  a  renormalization  similar  to  that  shown  in  Eq. 
(25). 
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Including  the  inhomogeneous  part  of  the  solution  with  the 
same  renormalization  as  shown  in  Eq.  (25),  we  find  that 
the  field  mode  amplitude  coefficients  at  the  front  interface 
of  the  grating  structure  are  given  by 


Qi(  0)  = 


'  c2a.,i  exp(-X2 a,,id) 

rd 

exp (-X2w>iz')X  Vi2t,ijYj(z')fe' , 

Jo  j 


'2o., i  > 


0  =s  i  <  2 p, 

2p  *£  i  <  4 p, 

(27) 


and  at  the  back  interface  of  the  grating  structure 
(z  —  d)  we  find  that 


'  c2w,i,  0  *£  i  <  2p, 

C2m,i  exp(k2o),irf)  +  exp(k2[u>id) 

X  f  exp(-k2W;i2')2  W2-i,^(z')d2', 

Jo  1 

2p  =s  i  <  4p. 

(28) 

Substituting  Eqs.  (27)  and  (28)  into  Eq.  (20),  we  can  ob¬ 
tain  the  space-harmonic  expansion  coefficients  of  the 
fields  at  the  second-harmonic  frequency  at  the  front  and 
back  interfaces  of  the  nonlinear  grating  layer. 


D.  Boundary  Condition  Matching 

To  obtain  the  complete  solution  of  the  fields  inside  the 
nanostructure,  it  is  necessary  to  determine  the  unknown 
coefficients  {c2w}  by  matching  the  tangential  components 
of  the  fields  at  each  interface.  The  fields  in  the  incident 
and  transmitted  half-spaces  on  either  side  of  the  grating 
can  be  represented  by  a  set  of  plane  waves  with  trans¬ 
verse  wave  vectors  matched  to  the  space-harmonic  modes 
inside  the  grating.  Thus  there  are  p  =  2a  +  1  plane 
waves  in  each  half-space,  with  complex  amplitude  coeffi¬ 
cients  collected  into  vectors  f  and  t  in  the  incident  and 
transmitted  half-spaces,  respectively.  The  space- 
harmonic  representation  of  the  transmitted  fields  at  the 
interface  ( z  =  d  1 )  can  be  written  as 


IS2lo,x(d+)\ 

S2„,y(d  +  ) 

U2„,x(d+) 

\*Wd+)/ 

where  I  is  the  identity  matrix  and  matrix  H  yields  the 
magnetic  field  amplitude  coefficients  from  the  electric- 
field  amplitude  coefficients  of  the  free-space  modes  in  the 
half-spaces  on  either  side  of  the  structure.  These  rela¬ 
tions  can  be  found  in  a  straightforward  way  by  using  Eq. 
(7),  and  an  explicit  definition  of  the  elements  of  matrix  H 
can  be  found  in  Eq.  (22)  or  Eq.  (26)  of  Ref.  20.  Similarly, 
the  space-harmonic  representation  of  the  fields  at  the  in¬ 
cident  boundary  ( z  =  0  )  is  given  by 
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where  the  vector  b  represents  the  incident  field.  For  the 
fundamental-frequency  waves,  the  vector  b  is  unity  for 
the  elements  corresponding  to  zeroth-order  electric  and 
magnetic  fields  and  zero  otherwise  [see  Eqs.  (21)  and  (22) 
of  Ref.  20] .  In  contrast,  for  the  second-harmonic  waves, 
the  field  is  assumed  to  be  generated  inside  the  nonlinear 
material,  and  the  externally  incident  field  is  uniformly 
zero. 

Owing  to  the  orthogonality  of  the  space-harmonic  basis, 
at  each  interface  every  transverse  space-harmonic  compo¬ 
nent  must  independently  satisfy  the  boundary  conditions. 
By  combining  Eqs.  (20),  (27),  and  (30)  for  the  fields  at  the 
incident  boundary,  we  obtain  the  expression 
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where  Q(ih)’  (0)  indicates  only  the  top  half  (0  =s  i 
<  2 p)  of  the  vector  Q(lh)(0),  and  the  expansion 
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has  been  applied.  Similarly,  by  combining  Eqs.  (20),  (28), 
and  (29)  for  the  fields  at  the  transmitted  boundary,  we  ob¬ 
tain  the  expression 
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where  Q(lh)’  ‘  (d)  indicates  only  the  bottom  half  (2 p  =S  i 
<  4 p)  of  the  vector  Q(lh)(d). 

Owing  to  numerical  stability  issues,  Eqs.  (31)  and  (32) 
cannot  be  solved  by  substitution.21  Instead,  we  combine 
them  into  a  single  system  of  equations: 
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Solution  of  Eq.  (33)  yields  the  complex  amplitude  coeffi¬ 
cients  r,  c2w,  and  f,  which  respectively  describe  the  re- 
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fleeted,  internal,  and  transmitted  fields  of  the  grating  at 
the  second-harmonic  frequency.  Thus  this  technique  per¬ 
mits  the  analysis  of  second-order  nonlinear  optical  phe¬ 
nomena  in  periodic  nanostructures  in  the  undepleted- 
pump  approximation. 


3.  NANOSTRUCTURE  DESIGN  FOR 
ENHANCED  SECOND-HARMONIC 
GENERATION 

A.  Transverse  Near-Field  Localization 

In  a  previous  paper8  we  described  the  transverse  localiza¬ 
tion  of  the  fundamental  frequency  field  in  a  periodic  sub¬ 
wavelength  nanostructure  and  performed  an  approximate 
analysis  of  the  enhancement  of  the  SHG  output  that  is 
due  to  this  localization.  This  approximate  analysis  was 
based  solely  on  transverse  effects,  and  thus  it  did  not  take 
into  account  the  mismatch  in  the  phase  velocities  of  the 
fundamental  and  second-harmonic  fields  or  the  effects  of 
the  finite  length  of  the  nanostructure.  For  a  more  accu¬ 
rate  analysis,  we  apply  the  extended  RCWA  method  de¬ 
scribed  above  to  investigate  a  similar  structure,  but  for  a 
very  small  grating  thickness  and  replacing  ultrashort- 
pulse  illumination  with  a  monochromatic  field.  The 
nanostructure  is  a  subwavelength  square  grating  (see 
Fig.  1)  of  period  A  =  0.65X,  fill  factor  F  =  0.09,  and 
depth  varying  from  d  —  0.005X  to  d  =  5. OX  in  units  of 
the  fundamental  wavelength  in  vacuum.  The  refractive 
index  of  the  material  is  taken  to  be  =  3.346  at  the  fun¬ 
damental  frequency  and  n2„,  =  3.539  at  the  second- 
harmonic  frequency,  corresponding  to  the  material  prop¬ 
erties  of  GaAs  at  1.907  pm  and  0.954  pm,  respectively. 
These  particular  values  were  selected  to  ensure  that  the 
absorption  coefficient  of  GaAs  is  negligible  for  both  the 
fundamental  and  the  second-harmonic  wavelengths. 
The  nonlinear  coefficient  is  taken  to  be  0.001  in  inverse 
units  of  the  field  amplitude,  corresponding  to  the  second- 
order  nonlinear  coefficient  of  GaAs  (y<2)  =  240  pm/V  from 


x 


Fig.  1.  Schematic  diagram  of  the  transverse  field  localization 
nanostructure  for  SHG  enhancement.  The  nanostructure  is  a 
subwavelength  periodic  square  grating  with  period  A  =  0.65X, 
fill  factor  F  =  0.09,  index  of  refraction  n„  =  3.346  at  the  funda¬ 
mental  frequency  and  n2al  =  3.539  at  the  second-harmonic  fre¬ 
quency,  and  depth  ranging  from  d  =  0.005X  to  d  =  5.00X  in 
units  of  the  fundamental  wavelength  in  vacuum.  The  structure 
is  assumed  to  be  infinite  and  periodic  in  the  x  direction  and  infi¬ 
nite  and  uniform  in  the  y  direction.  The  nanostructure  is  illu¬ 
minated  by  a  normally  incident  plane  wave  at  the  fundamental 
frequency  with  wave  vector  kinc  and  electric  field  Einc  polarized  in 
the  y  direction. 


Fig.  2.  Total  transmitted  SHG  output  intensity  for  the  nano¬ 
structure  shown  in  Fig.  1  and  bulk  nonlinear  material  as  a  func¬ 
tion  of  depth. 


Ref.  15)  with  an  applied  field  amplitude  of  4.17 
X  106  V/m.  The  fundamental  and  second-harmonic 
fields  are  both  assumed  to  be  polarized  in  the  y  direction, 
as  shown  in  Fig.  1.  A  nanostructure  of  this  type  would  be 
fabricated  most  easily  either  by  deep  etching  into  a  sub¬ 
strate  or  a  multilayer  growth;  the  resulting  structure 
would  differ  from  that  shown  in  Fig.  1  because  of  the  pres¬ 
ence  of  the  substrate  and  the  finite  extent  of  the  structure 
in  the  transverse  direction.  Although  it  is  possible  to 
perform  the  analysis  for  a  structure  that  incorporates 
these  characteristics,  the  interpretation  of  the  results 
would  be  significantly  more  complicated  and  therefore 
less  clear.  Thus  we  analyze  the  simple  nanostructure 
shown  in  Fig.  1 — which  is  assumed  to  be  infinitely  peri¬ 
odic  in  the  x  direction,  infinite  and  uniform  in  the  y  direc¬ 
tion,  finite  and  uniform  in  the  z  direction,  surrounded  by 
vacuum  on  both  sides,  and  not  explicitly  including  the 
substrate — in  order  to  investigate  and  present  in  a 
straightforward  manner  the  fundamental  characteristics 
of  the  optical  field  interaction  with  the  nanostructure. 

Using  the  RCWA-based  undepleted-pump  SHG  model¬ 
ing  tool,  we  investigate  the  total  transmitted  SHG  output 
of  the  nanostructure  and  compare  it  with  that  of  the  bulk 
nonlinear  material.  Figure  2  shows  the  total  transmit¬ 
ted  SHG  output  of  the  structure  as  a  function  of  depth  for 
the  nanostructure  and  the  bulk  material.  For  a  certain 
range  of  depths,  the  nanostructure  yields  a  greater  SHG 
intensity  output  than  the  bulk  material.  The  peak  oc¬ 
curs  for  a  depth  of  d  =  0.29X,  with  an  output  from  the 
nanostructure  approximately  2.05  times  that  of  the  bulk 
material  of  the  same  thickness.  For  greater  material 
thicknesses,  however,  in  most  cases  the  nonlinear  output 
of  the  bulk  material  exceeds  that  of  the  nanostructure. 
The  oscillations  in  the  total  transmitted  SHG  intensity  as 
a  function  of  depth  for  both  the  nanostructure  and  the 
bulk  material  are  due  to  Fabry-Perot  effects  and  phase 
velocity  mismatch  between  the  fundamental  and  second- 
harmonic  waves  and  will  be  considered  in  the  following. 

In  order  to  analyze  further  the  nonlinear  output  of 
the  nanostructure,  it  is  useful  first  to  consider  the  results 
for  the  bulk  nonlinear  material  shown  in  Fig.  2.  It  is 
clear  that  the  transmitted  SHG  output  as  a  function  of 
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the  nonlinear  material  thickness  results  from  the  super¬ 
position  of  several  phenomena.  First,  the  thin  slab  of 
nonlinear  material  acts  as  a  Fabry-Perot  cavity  owing  to 
the  reflection  from  the  front  and  back  interfaces,  thus 
modulating  the  output  of  the  SHG  process.  For  a  given 
vacuum  wavelength  X',  the  Fabry-Perot  resonances  are 
found  at  integer  multiples  of  l  =  \'/(2 n).  The  lowest- 
order  resonant  material  thicknesses  for  the  fundamental 
and  second-harmonic  frequencies  are  /,„  =  0.149X  and 
=  0.071X,  respectively.  Although  an  approximation, 
it  can  be  seen  in  Fig.  2  that  the  integer  multiples  of  these 
two  characteristic  lengths  correspond  to  the  positions  of 
the  numerous  narrow  peaks  in  the  total  transmitted  SHG 
of  the  bulk  material.  Second,  we  must  also  consider  the 
phase  mismatch  between  the  fundamental  and  second- 
harmonic  fields  that  is  due  to  material  dispersion.  The 
effective  propagation  speeds  of  the  fields  inside  the  slab  at 
the  fundamental  and  second-harmonic  frequencies  are 
v,„  =  0.2989c  and  v2w  —  0.2826c,  respectively  (where  c  is 
the  speed  of  light  in  vacuum),  corresponding  to  a  charac¬ 
teristic  length  of  lc  =  1.30X.  Since  this  is  the  distance 
for  the  fundamental  and  second-harmonic  fields  to 
dephase  by  tt,  we  expect  the  SHG  output  to  reach  a  maxi¬ 
mum  at  odd  multiples  of  this  characteristic  length  and  to 
return  to  zero  at  even  multiples  of  this  length  as  a  result 
of  destructive  interference.  This  corresponds  to  the  en¬ 
velope  of  the  total  transmitted  SH  intensity  function  for 
the  bulk  material  case  as  seen  in  Fig.  2. 

To  understand  the  SHG  output  characteristics  of  the 
nanostructure,  we  will  perform  a  similar  analysis.  How¬ 
ever,  we  must  first  determine  the  propagation  character¬ 
istics  of  the  nanostructure  modes.  Using  the  RCWA 
method,  we  can  obtain  this  information  from  the  eigenval¬ 
ues  and  eigenvectors  of  the  periodic  structure  [see  Eq. 
(A8)].  Each  eigenvector  describes  one  of  the  natural 
modes  of  the  nanostructure,  and  the  corresponding  eigen¬ 
value  describes  its  propagation  characteristics — the  real 
part  represents  the  decay  coefficient,  and  the  imaginary 
part  represents  the  propagation  constant  (which  is  analo¬ 
gous  to  the  wave  vector  for  a  plane  wave).  Furthermore, 
owing  to  the  symmetry  of  the  nanostructure  with  respect 
to  the  z  axis  (the  finite  dimension  of  the  nanostructure), 
the  eigenmodes  will  exist  in  pairs  with  eigenvalues  of  op¬ 
posite  sign,  representing  equivalent  modes  propagating  in 
opposite  directions.  For  clarity,  in  the  following  only  one 
of  each  pair  of  eigenmodes  is  shown.  Figure  3  shows  re¬ 
sults  of  an  analysis  of  the  nonevanescent  energy-carrying 
modes  of  the  nanostructure,  obtained  by  selecting  only 
those  eigenmodes  with  a  real  part  of  the  eigenvalue  less 
than  10  s  in  magnitude  and  an  amplitude  coefficient 
greater  than  0.01  in  magnitude.  At  the  fundamental 
frequency  there  is  only  a  single  mode,  shown  in  Fig. 
3(a),  which  is  transversely  localized  in  the  high- 
refractive-index  region  of  the  nanostructure,  as  predicted 
in  our  previous  work.8  In  addition,  there  are  two  modes 
at  the  second-harmonic  frequency,  shown  in  Figs.  3(b)  and 
3(c),  one  localized  in  the  high-refractive-index  region  and 
one  localized  primarily  in  the  air  gap  of  the  nanostruc¬ 
ture.  Figure  3  also  shows  the  real  and  imaginary  parts 
of  the  eigenvalues  for  each  of  the  modes.  These  con¬ 
stants  correspond  to  Fabry-Perot  resonance  characteris¬ 
tic  lengths  of  lw  =  0.269X  for  the  fundamental  mode  and 
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Fig.  3.  Properties  of  the  nonevanescent  energy-carrying  eigen¬ 
modes  of  the  nanostructure  shown  in  Fig.  1  for  the  (a)  fundamen¬ 
tal  and  (b),  (c)  second-harmonic  wavelengths  computed  by  using 
the  modified  RCWA  tool:  For  each  mode  the  transverse  profile 
in  one  period  of  the  nanostructure  as  well  as  the  attenuation  co¬ 
efficient  and  propagation  constant  (real  and  imaginary  parts  of 
the  corresponding  eigenvalue  Z  respectively)  are  shown. 


~  0.278X  and  l2lo2  =  0.092X  for  the  second- 
harmonic  modes.  The  phase  velocities  of  the  three  modes 
shown  in  Figs.  3(a)— 3(c)  are  v m  =  0.538c,  v2w  i  —  1.10c, 
and  u2<u,2  =  0.370c,  corresponding  to  phase  matching 
characteristic  lengths  of  lc  l  =  0.527X  and  lc2  =  0.592X. 
These  parameters  describe  qualitatively  the  transmitted 
SHG  intensity  as  a  function  of  depth  from  the  nanostruc¬ 
ture. 

For  this  nanostructure  geometry,  phase  matching  be¬ 
tween  the  fundamental  and  either  of  the  second-harmonic 
modes  is  not  achieved  for  these  depth  ranges.  Further¬ 
more,  since  the  difference  in  phase  velocity  is  greater  for 
the  nanostructure  than  for  the  bulk  material,  we  do  not 
expect  in  general  to  be  able  to  achieve  enhanced  SHG  out¬ 
put  from  this  nanostructure  relative  to  the  bulk  material, 
despite  the  transverse  field  localization.  In  fact,  the  par¬ 
ticular  depths  discussed  above  where  the  nanostructure 
yields  a  greater  SHG  transmitted  intensity  than  the  bulk 
material  are  due  to  the  differing  effective  indices  of  re¬ 
fraction  and  the  Fabry-Perot  modes  of  the  structure,  re¬ 
sulting  in  a  significantly  higher  fundamental  field  ampli¬ 
tude  in  the  nanostructure  in  comparison  with  the  bulk 
material  of  the  same  depth.  However,  the  peak  transmit¬ 
ted  SH  intensity  from  the  nanostructure  is  approximately 
30%  of  that  obtained  from  the  bulk  material,  although  the 
nanostructure  has  only  9%  of  the  nonlinear  material  by 
volume  as  well  as  a  greater  phase  mismatch  between  the 
fundamental  and  SH  fields.  This  fact  suggests  that  the 
nanostructure  does  in  fact  contribute  to  the  enhancement 
of  the  SHG  process  through  transverse  localization  of  the 
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field;  but  to  achieve  superior  performance,  phase  match¬ 
ing  must  be  realized  in  the  nanostructure. 

B.  Phase  Matching 

To  implement  a  SHG  enhancement  nanostructure  that  in¬ 
corporates  both  transverse  field  localization  and  phase 
matching,  we  investigate  an  alternative  nanostructure 
having  two  ridges  per  period,  instead  of  one.  In  Subsec¬ 
tion  3. A  we  found  that  the  fundamental  frequency  mode 
of  the  nanostructure  is  transversely  localized  in  the  high- 
refractive-index  region,  while  the  second-harmonic  fre¬ 
quency  mode  of  interest  is  transversely  localized  in  the 
low-refractive-index  region  and  has  a  phase  velocity  sig¬ 
nificantly  higher  than  that  of  the  fundamental  frequency 
mode.  Thus  we  expect  that  the  introduction  of  a  second, 
smaller  nonlinear  material  ridge  into  the  center  of  the 
low-refractive-index  region  will  affect  the  phase  velocity 
of  the  second-harmonic  mode  more  dramatically.  Specifi¬ 
cally,  the  introduction  of  a  high-refractive-index  material 
in  a  region  where  the  fundamental  mode  has  a  low  ampli¬ 
tude  but  the  second-harmonic  mode  has  a  large  amplitude 
will  serve  to  significantly  reduce  the  phase  velocity  of  the 
second  harmonic  mode  while  leaving  the  fundamental 
mode  largely  unaffected.  By  engineering  this  differential 
slowing  of  the  two  modes  to  exactly  compensate  for  the 
existing  phase  velocity  mismatch  between  the  two  modes, 
we  can  achieve  phase  matching.  It  is  important  to  point 
out  that  the  second  material  ridge  must  be  smaller  than 
the  original  nanostructure  ridge;  otherwise,  the  mode 
properties  of  the  nanostructure  would  be  qualitatively 
changed.  As  in  the  previous  example,  for  clarity  we  con¬ 
sider  a  nanostructure  that  is  infinitely  periodic  in  the  x 
direction,  infinite  and  uniform  in  the  y  direction,  finite, 
uniform  in  the  z -direction,  surrounded  by  vacuum  on  both 
sides,  and  not  explicitly  including  the  substrate. 

A  schematic  drawing  of  the  structure  is  shown  in  Fig.  4. 
The  period  of  the  grating  is  again  A  =  0.65X,  and  the  fill 
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Fig.  4.  Schematic  diagram  of  the  modified  transverse  field  local¬ 
ization  nanostructure  with  improved  phase  matching,  consisting 
of  a  periodic  subwavelength  grating  having  two  nonlinear  mate¬ 
rial  ridges  with  fill  factors  F1  =  0.063  and  F2  =  0.03  per  period 
A  =  0.65X.  The  index  of  refraction  is  n,  =  3.346  at  the  funda¬ 
mental  frequency  and  n2al  =  3.539  at  the  second-harmonic  fre¬ 
quency,  and  the  depth  ranges  from  d  =  0.005X  to  d  =  5.00X  in 
units  of  the  fundamental  wavelength  in  vacuum.  The  structure 
is  assumed  to  be  infinite  and  periodic  in  the  x  direction  and  infi¬ 
nite  and  uniform  in  the  y  direction.  The  nanostructure  is  illu¬ 
minated  by  a  normally  incident  plane  wave  at  the  fundamental 
frequency  with  wave  vector  kinc  and  electric  field  Einc  polarized  in 
the  y  direction. 
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Fig.  5.  Properties  of  the  nonevanescent  energy-carrying  eigen- 
modes  of  the  modified  nanostructure  shown  in  Fig.  4  for  the  (a), 
(b)  fundamental  and  (c),  (d)  second-harmonic  wavelengths  com¬ 
puted  by  using  the  modified  RCWA  tool:  For  each  mode  the 
transverse  profile  in  one  period  of  the  nanostructure  as  well  as 
the  attenuation  coefficient  and  propagation  constant  (real  and 
imaginary  parts  of  the  corresponding  eigenvalue  Z)  are  shown. 


factors  of  the  two  ridges  are  F1  =  0.063  and  F2  =  0.03. 
Analyzing  the  nonevanescent  energy-carrying  modes  of 
the  nanostructure,  we  find  two  modes  each  at  the  funda¬ 
mental  and  second-harmonic  frequencies,  which  are 
shown  in  Figs.  5(a)— 5(d).  At  each  frequency,  one  mode  is 
localized  in  the  wide  nonlinear  material  ridge  (with  fill 
factor  F x)  shown  in  the  center  of  the  figure,  while  the 
other  mode  is  localized  primarily  in  the  narrow  second 
ridge  (with  fill  factor  F 2).  Figure  5  also  shows  the  real 
and  imaginary  parts  of  the  mode  propagation  constants  of 
the  four  modes. 

To  verify  the  ability  to  control  the  phase  velocities  of 
the  fundamental  and  second-harmonic  modes  by  using 
the  geometry  of  the  nanostructure,  we  investigate  the  de¬ 
pendence  of  the  mode  phase  velocities  on  the  fill  factors 
F x  and  F2  of  the  nanostructure.  We  compute  the  phase 
velocities  of  the  four  modes  as  a  function  of  the  narrow- 
ridge  fill  factor  F 2  while  keeping  the  wide-ridge  fill  factor 
F1  =  0.063  constant,  as  shown  in  Fig.  6.  In  this  case,  the 
fundamental  and  second-harmonic  modes  that  are  local¬ 
ized  in  the  wide  F1  ridge  [labeled  F-l  and  SH-1  in  Fig.  6, 
corresponding  to  the  modes  shown  in  Figs.  5(a)  and  5(c), 


1926  J.  Opt.  Soc.  Am.  A/Vol.  19,  No.  9/September  2002 


0.01  0.02  0.03  0.04  0.05 

Second  Ridge  Fill  Factor  F2 


Fig.  6.  Phase  velocities  of  the  two  fundamental  (F-l  and  F-2) 
and  two  second-harmonic  (SH-1  and  SH-2)  nonevanescent 
energy-carrying  eigenmodes  of  the  nanostructure  shown  in  Fig.  4 
as  a  function  of  fill  factor  F2  with  F1  =  0.063  constant. 


Fig.  7.  Total  transmitted  SHG  output  intensity  for  the  modified 
nanostructure  shown  in  Fig.  4  and  bulk  nonlinear  material  as  a 
function  of  depth. 


respectively]  have  an  approximately  constant  phase  ve¬ 
locity  with  respect  to  F2 ,  while  the  other  two  modes  [la¬ 
beled  F-2  and  SH-2  in  Fig.  6,  corresponding  to  the  modes 
shown  in  Figs.  5(b)  and  5(d),  respectively]  exhibit  a  de¬ 
pendence  of  their  phase  velocity  on  the  parameter  F2  ■ 
Thus  the  phase  velocities  of  the  fundamental  mode  local¬ 
ized  in  the  wide  F 1  ridge  and  the  second-harmonic  mode 
localized  in  the  narrow  F 2  ridge  can  be  adjusted  essen¬ 
tially  independently  through  the  geometry  of  the  nano¬ 
structure.  More  important,  the  phase  velocities  can  be 
made  equal  to  achieve  phase  matching  for  SHG,  as  shown 
in  Fig.  6,  for  the  nanostructure  geometry  shown  in  Fig.  4. 

Using  the  mode  propagation  constants  from  Figs.  5(a), 
5(b),  5(c),  and  5(d)  for  the  geometry  shown  in  Fig.  4,  we 
find  the  phase  velocities  of  the  modes  to  be  uml 
=  0.630c,  0^,2  =  1.32c,  u2w,i  =  0.428c,  and  v2l0,2 

=  0.630c,  respectively.  The  addition  of  the  second 
high-refractive-index  ridge  results  in  good  matching  of 
the  phase  velocities  of  the  “slow”  fundamental  frequency 
mode  (o^i)  and  the  “fast”  second-harmonic  frequency 
mode  (o2W;2).  The  characteristic  length  of  the  propaga¬ 
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tion  constant  mismatch  of  these  two  modes  is  lc 
—  981.84.  As  a  result,  we  expect  that  phase  mismatch 
will  not  significantly  affect  the  output  of  the  nonlinear 
process  in  the  modified  nanostructure,  at  least  over  the 
practical  nanostructure  sizes  under  consideration. 
Phase  matching  in  the  modified  nanostructure  results  in 
a  significant  enhancement  of  the  second-harmonic  output 
relative  to  the  original  nanostructure  design  as  well  as  to 
the  bulk  nonlinear  material. 

A  comparison  of  the  total  transmitted  SHG  intensity 
of  the  modified  nanostructure  and  bulk  material  as 
a  function  of  the  nonlinear  material  depth  is  shown  in 
Fig.  7.  For  both  the  nanostructure  and  the  bulk  mate¬ 
rial,  Fabry-Perot  effects  still  strongly  affect  the  SHG  out¬ 
put.  However,  maxima  of  the  total  transmitted  SHG  out¬ 
put  of  the  nanostructure  are  significantly  higher  than 
those  of  the  bulk  material.  For  example,  for  a  depth 
d  =  0.3954,  the  nanostructure  yields  a  total  transmitted 
SHG  intensity  of  approximately  2.4  times  the  maximum 
obtained  from  the  phase-mismatched  bulk  material  for 
any  depth  shown;  and  for  a  depth  d  =  2.3824,  the  ratio  is 
approximately  463.  As  shown  in  Fig.  7,  the  peak  corre¬ 
sponding  to  the  enhanced  SHG  output  in  the  latter  case  is 
extremely  narrow,  suggesting  that  the  elevated  output  is 
strongly  dependent  on  the  longitudinal  resonances  of  the 
nanostructure.  In  contrast,  in  the  former  case  the  width 
of  the  enhanced  SHG  output  peak  is  relatively  broad,  sug¬ 
gesting  that  the  lower  enhancement  factor  can  be 
achieved  without  extreme  sensitivity  to  the  nanostruc¬ 
ture  geometry  (thus  requiring  less  rigid  fabrication  toler¬ 
ances).  In  addition,  in  the  former  case  the  structure 
depth  is  significantly  less  than  the  characteristic  phase 
matching  length  in  the  bulk  material,  indicating  that  the 
elevated  SHG  output  of  the  modified  nanostructure  is  not 
solely  the  result  of  longitudinal  resonances  or  improved 
phase  matching  in  the  nanostructure.  Furthermore,  the 
fact  that  the  modified  nanostructure  consists  of  only  9.3% 
nonlinear  material  by  volume  emphasizes  the  importance 
of  the  transverse  near-field  localization  effect  in  the  nano¬ 
structure.  Thus  these  results  demonstrate  that  sub¬ 
wavelength  periodic  optical  nanostructures  can  be  ap¬ 
plied  to  enhance  SHG  via  two  mechanisms:  (1) 
transverse  field  localization  in  the  nonlinear  material  to 
increase  the  peak  power  and  (2)  engineering  of  the  propa¬ 
gation  speeds  of  the  fundamental  and  SH  modes  to 
achieve  phase  matching. 

4.  CONCLUSION 

We  have  extended  the  well-established  RCWA  algorithm 
to  analyze  SHG  in  periodic  nanostructures  in  the 
undepleted-pump  approximation.  This  nonlinear  elec¬ 
tromagnetic  modeling  tool  permits  rigorous  analysis  of 
the  structure  and  propagation  characteristics  of  the  nano¬ 
structure  modes,  facilitating  the  investigation  of  near¬ 
field  phenomena  in  the  nanostructures  as  well  as  the  de¬ 
sign  of  novel  nonlinear  optical  devices.  In  principle,  it 
would  be  possible  to  continue  the  perturbation  expansion 
of  the  nonlinear  polarization  term,  permitting  the  analy¬ 
sis  of  downconversion  (including  the  depleted-pump  case) 
as  well  as  other  nonlinear  optical  phenomena  of  interest, 
but  much  more  computational  power  would  be  required. 
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In  addition,  using  this  tool,  we  have  investigated  two 
nanostructures  for  enhanced  SHG.  In  the  first  case,  a 
structure  for  transverse  localization  of  the  field  in  the 
nonlinear  material  provides  evidence  of  enhanced  perfor¬ 
mance  resulting  from  the  field  localization,  but  the  overall 
transmitted  SHG  output  of  the  device  was  lower  than 
that  of  the  bulk  material  owing  to  phase  mismatch  in  the 
structure.  In  the  second  case,  the  nanostructure  was 
modified  to  achieve  phase  matching  through  the  introduc¬ 
tion  of  a  second  high-refractive-index  material  ridge  per 
period.  This  structure  yields  increased  SHG  transmitted 
intensity  in  comparison  with  the  bulk  nonlinear  material, 
more  than  two  orders  of  magnitude  in  some  cases,  despite 
the  fact  that  it  is  composed  of  only  9.3%  nonlinear  mate¬ 
rial  by  volume.  Thus,  using  the  additional  design  de¬ 
grees  of  freedom  provided  by  the  periodic  nanostructure 
makes  it  possible  to  construct  a  device  capable  of  enhanc¬ 
ing  SHG  efficiency  through  both  transverse  near-held  lo¬ 
calization  and  phase  matching  in  an  isotropic  nonlinear 
material.  As  a  result,  this  approach  could  provide  sig¬ 
nificantly  improved  performance  in  comparison  with 
other  microstructure-  and  nanostructure-based  tech¬ 
niques  to  enhance  SHG  efficiency,  albeit  by  using  a  more 
complex  nanostructure  with  smaller  features.  The  rela¬ 
tive  merits  of  the  various  optical  nonlinearity  enhance¬ 
ment  techniques  depend  strongly  on  the  nature  and  re¬ 
quirements  of  the  specific  application  of  interest;  thus  a 
direct  comparative  analysis  is  beyond  the  scope  of  this  pa¬ 
per  and  will  be  investigated  in  the  future. 

Finally,  the  presence  of  significant  Fabry-Perot  effects 
indicates  that  longitudinal  localization  of  the  held  in  the 
cavity  formed  by  the  nanostructure  is  present  and  dra¬ 
matically  inhuences  the  SHG  output.  In  fact,  the  appli¬ 
cation  of  the  longitudinal  modes  of  resonant  cavities  to 
enhance  the  held  amplitude  and  thus  nonlinear  optical  ef¬ 
fects  is  a  well-known  technique.  More-advanced  device 
designs  could  eliminate  cavity  effects  through  nanostruc- 
tured  antireflection  layers  to  relax  design  tolerances  or  to 
take  advantage  of  cavity  phenomena  by  using  photonic 
crystal  techniques  to  further  enhance  the  nonlinear  out¬ 
put.  Finally,  an  important  beneht  of  this  approach  to  en¬ 
hance  the  efficiency  of  SHG  in  nanostructures  is  that  the 
materials  and  designs  used  are  compatible  with  standard 
microfabrication  techniques,  facilitating  integration  with 
other  photonic  devices  and  systems. 


APPENDIX  A:  SUMMARY  OF  STANDARD 
RIGOROUS  COUPLED-WAVE  ANALYSIS 
FORMULATION 

This  appendix,  included  for  completeness,  summarizes 
the  standard  formulation  of  RCWA  and  introduces  the  no¬ 
tation  used  in  this  paper.  A  more  detailed  description  of 
the  method  can  be  found  in  Refs.  20-24.  The  grating  re¬ 
gion  is  taken  to  be  infinitely  periodic  in  one  transverse  di¬ 
mension  ( x )  with  grating  period  A,  infinite  and  uniform  in 
the  other  transverse  dimension  (y),  and  finite  and  uni¬ 
form  in  the  depth  dimension  (z).  The  dielectric  constant 
function  inside  the  grating  may  be  a  tensor  to  account  for 
anisotropic  materials,  but  it  depends  only  on  the  coordi¬ 
nate  x  and  is  a  periodic  function: 


e(r)  =  e(r  +  Ax) 


|  ejx),  0  <  z  <  d 
[  e0  z  <  0;z  >  d 


(Al) 


The  incident  field  is  assumed  to  be  a  monochromatic 
plane  wave  of  unit  amplitude  with  wave  vector  kinc .  Un¬ 
der  these  assumptions,  the  electric  and  magnetic  fields  in¬ 
side  the  grating  region  may  be  represented  by  using  a 
space-harmonic  field  expansion: 


Ea,(r)  =  2  {Su,i,x (z)x  +  Sa'iiy(z)$ 

i 

+  SWtitZ(z)z}exTp(-jd(0ti  ■  r),  (A2) 

HJr,  t)  =  A/—  2  [U^i^x  +  U^i  (z)y 
v  Mo  i 

+  UWji,z(z)2]exp(-y'd:IU;i  •  r),  (A3) 

where  Swi  and  are  the  space-harmonic-expansion 
amplitude  coefficients  for  the  electric  and  magnetic  fields, 
respectively,  e0  is  the  dielectric  permittivity  of  free  space, 
and  /in  is  the  magnetic  permeability  of  free  space.  The 
space-harmonic  expansion  contains  a  total  of  p  =  2a 
+  1  elements,  referenced  by  index  i  e  {— a,...,  —  1, 
0,  l,...,a}.  The  transverse  component  of  the  wave  vec¬ 
tor  is  given  by 

i  i^inc,x  lKg)x  "1  &inc,yy >  (A4) 

where  Kg  =  2  tt/A  .  Since  the  grating  is  periodic,  the 
{vw)  component  ( v ,  w  <e  {x,  y,  z})  of  the  dielectric  con¬ 
stant  tensor  can  also  be  expanded  in  a  space-harmonic  se¬ 
ries: 


em,vw(x)  =  eo2  G,i,»«.exp (jkKgx).  (A5) 

k 


A  similar  expansion  is  performed  for  the  electric  and  mag¬ 
netic  fields  in  the  half-spaces  on  either  side  of  the  grating 
region.  Inserting  Eqs.  (A2),  (A3),  and  (A5)  into  Maxwell’s 
equations  and  simplifying  produces  a  coupled  system  of 
first-order  partial  differential  equations: 


S«vc(2)  \ 

jsa 

*(z) 

S^Z) 

=  Aw 

^  (o,y(z  ) 

uw,x(z )  1 

Az> 

U^y(z)/ 

\  UWty(z)l 

where  the  matrix  AnJ  is  determined  by  the  structure  of  the 
grating  and  Maxwell’s  equations  [see,  for  example,  Eq. 
(32)  of  Ref.  24  for  a  detailed  description],  and  the  vector 
Sv(z)  contains  the  p  =  2a  +  1  elements  of  the  space- 
harmonic  expansion  coefficients 


^G),*(z) 


Sc,,-a(z)  \ 

^  (o,—  a  +  l(^  ) 

Sa,a(z)  I 


(A7) 


with  a  similar  expansion  for  the  vectors  S^^iz),  U„l>x(z), 
and  U„jy(z). 

Using  the  eigenvector  matrix  W„,  corresponding  to  the 
matrix  A„, ,  we  can  write  the  solution  to  Eq.  (A6)  as 
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S»,x(z)\ 
^  w,y  ) 

Um,x{z) 

\u«,,y(z)l 


W,„  expi'Dz  )cm , 


(A8) 


where  cw  is  a  vector  of  length  4 p  containing  the  ampli¬ 
tude  coefficients  of  each  eigenmode  and  the  matrix  D  is  a 
diagonal  matrix  containing  the  eigenvalues.  Using  Eq. 
(A8)  to  propagate  the  fields  across  the  grating  region,  and 
applying  the  appropriate  boundary  conditions  at  each  in¬ 
terface,  we  can  find  the  amplitude  coefficients  cw  (see,  for 
example,  Ref.  20  for  a  more  detailed  explanation). 
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